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Abstract 

The couplings of pions with heavy baryons 52(^*7 S) and <?3(£*,A) are studied 
with light-cone QCD sum rules in the leading order of heavy quark effective theory. 
Both sum rules are stable. Our results are #2 = 1.56±0.3±0.3, 53 = 0.94±0.06±0.2. 
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^ ! 1 Introduction 

a: 

Important progress has been achieved in the interpretation of heavy hadrons composed of a 
heavy quark with the development of the heavy quark effective theory (HQET) [|TJ . HQET 
provides a systematic expansion of the heavy hadron spectra and transition amplitude 
in terms of I/tuq, where tjiq is the heavy quark mass. Of course one has to employ 
some specific nonperturbative methods to arrive at the detailed predictions. Among the 
various nonpeturbative methods, QCD sum rules is useful to extract the low-lying hadron 
properties 0. 

The couplings of the heavy mesons with pions has been analysed with QCD sum rules 
[5-15]. The couplings of heavy baryons with soft pions are estimated from QCD sum rules 
in an external axial field [|lj|] . In this approach the mass difference A between the baryons 
in the initial and final states is approximately taken to be zero. 

In this work we employ the light-cone QCD sum rules (LCQSR) in HQET to calculate 
the couplings #2,3 to the leading order of l/m,Q. The LCQSR is quite different from the 
conventional QCD sum rules, which is based on the short-distance operator product ex- 
pansion (OPE). The LCQSR is based on the OPE on the light cone, which is the expansion 
over the twists of the operators. The main contribution comes from the lowest twist oper- 
ator. Matrix elements of nonlocal operators sandwiched between a hadronic state and the 
vacuum defines the hadron wave functions. The LCQSR approach has the advantage that 
the double Borel transformation is used so that the the continuum contribution is treated 
in a way better than the external field approach. Moreover, the final sum rule depends 
only on the value of the wave function at a specific point like ^p n (u = |), which is much 
better known than the whole wave function 
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2 Sum rules for the coupling constants 

We first introduce the interpolating currents for the heavy baryons: 

Va(x) = t abc [u aT {x)C lb d\x)}h c v {x) , (1) 
Vx+(x) = ^{u^ix^^ixM^Kix) , (2) 

= e abc [u aT (x)C lu u b (x)](- 9 r + \lH)K{*) , (3) 

where a, b, c is the color index, u(x), d(x), h v (x) is the up, down and heavy quark fields, T 
denotes the transpose, C is the charge conjugate matrix, g^ v = g^ v — v^v", jjf = 7 M — vv M , 
and i> M is the velocity of the heavy hadron. 

The overlap amplititudes of the interpolating currents with the heavy baryons is defined 

as: 

(0|tJ A |A)=/AUA, (4) 

(0|7fc|£) = / E tt E , (5) 
(0|^|S*) = ^,, (6) 
where w^* is the Rarita-Schwinger spinor in HQET. In the leading order of HQET, / E = / E * 
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We adopt the same notations for g 2 ,3 as in |I9| . The coupling constants gi and are 
defined through the following amplitudes: 

M(E: - A c tt) = ifu Ac q^ t , (7) 

JIT 

M(E* S c tt) = ^j^ie vpa X<lt^oltlbU^ c , (8) 

where = 132MeV, g M is the pion momentum. The process @ is kinematically forbidden. 
It is very important to get a reliable estimate of the coupling g 2 since it is not directly 
accessible experimentally. 

In order to derive the sum rules for the coupling constants we consider the correlators 

fd*x e-*-(7r( g )|r(i^.(0)^(x)) |0> = ^{-gT + \lH) (9) 

d 4 x e- ik -*{ir(q)\T(ri^(0)fj A (x)) |0) = 1±^^(-^ + ph^A","') , (10) 

where k' = k — q, q 1 ^ = q^ — (q ■ v)v^, uj = 2v ■ k, uj' = 2v ■ k! and q 2 = 0. 

Let us first consider the function Ge*a(^,u/) in (0). As a function of two variables, 
it has the following pole terms from double dispersion relation 

^ + T ' + / , (11) 



V^fn (2A S * - uj') (2A A - w) 2A S * - uj' 2A A - w 
where / E * etc are constants defined in A s * = m E » — mq. 
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Neglecting the four particle component of the pion wave function, the expression for 
G-2* t \(u),u') with the tensor structure reads 

4^°° dt J dxe~ ikx S(-x - vt)Tr{(7r(q)\u{0)d{x)\0)[^CiS T {-x)C lu }} , (12) 

where iS(—x) is the full light quark propagator with both perturbative term and contri- 
bution from vacuum fields. 



x 



iS(x) = (0\T[q(x),q(Om 
' " — — (qg s a ■ Gq) 



167T 2 JO 



2tt 2 x 4 12 192 
du{^-a ■ G{ux) - Aiu^rG^(ux) ln u] + 



(13) 



x z 



Similarly for G-%* p(ui , u') we have: 

4 jf dt J dxe~ ikx 5{-x - vt)Tr{{it{q)\u{Q)d{x)\Q)[ lp CiS T {-x)C lv }} , (14) 
To the present approximation, we need the following OPE near the light cone for two- 



and three-particle pion wave functions [12 



< n(q)\ rf»7 M 75«(0)|0 > = \ du e™*») + x 2 9l (u) + G(x 4 )) 

J 

2 n rl 

due iuqx g 2 (u) 
du e iuqx y P (u) , 



< 7r(g)|d(x)* 75 w(0)|0> 



qx Jo 



m u + m d Jo 



< 7r(q)\d(x)a fJ ,ul5u(0)\0 > = i(q^x u - q v x^ 



6(m u + m d ) Jo 



du e iuqx ^ a {u) 



(15) 
(16) 
(17) 



< n(q)\d(x)a a/3 j5g s G pu (ux)u(0)\0 >= 

if-in[{qf,q a gvf3 - q v q a 9tf) ~ {%%9v a - ^<?/?fiw)] / ^{a t )e iqx{ax+vaz) , ill 



and 



< 7r(g)|d(x)7 M 755f s G a/3 (t;a;)n(0)|0 > 



qp\9^ - —7) - q a \9fo 



q ■ x 



q ■ x 



Vai(p±(ai)e 



\ iqx(a 1 +va 3 ) 



+U—(qaXf3-qpx a ) I Vaiip^e^' "" ' ' 



< n(q)\d(x)'j IM g s G af3 (vx)u(0)\0 > 



fn[qp(g*n~^z) -q a (. 



q ■ x 



9fo 



q ■ x ' 



)1 f Da^ifa.le" 1 '^"" 3 ' 



q ■ x 



(19) 



(20) 
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The operator G a p is the dual of G a $: G a p = \e a pspG sp ; Van is defined as Voti = 
daida2da35(l—ai—a2—a3). Due to the choice of the gauge x^A^x) = 0, the path-ordered 
gauge factor Pexp (ig s Jq dux p A^{ux)) has been omitted. 

The wave function ip n (u) is associated with the leading twist 2 operator, gi(u) and 
g 2 (u) correspond to twist 4 operators, and (fp(u) and <p a {u) to twist 3 ones. The function 
y? 37r is of twist three, while all the wave functions appearing in eqs.([T9|), (p0|) are of twist 
four. The wave functions if(xi, fi) (fi is the renormalization point) describe the distribution 
in longitudinal momenta inside the pion, the parameters Xi (J2i x i — 1) representing the 
fractions of the longitudinal momentum carried by the quark, the antiquark and gluon. 

The wave function normalizations immediately follow from the definitions (JT5|)-([ZU|): 
fodu ip n (u) = fodu (p a (u) = 1, Jq du gi(u) = S 2 /12, jVa^^ai) = J Va^w (a») = 0, 
/ Va>i(p±(a>i) = — J Uai<p\\(ai) = 6 2 /3, with the parameter 6 defined by the matrix element: 
< n(q)\dg s G a ^ a u\0 >= i5 2 f ir q^ t . 



Expressing (|T2| ) and (|14|) with the pion wave functions, we arrive at: 

i . r°° . r 1 . .v, .Awt -u't r 1 



i poo rl . i 

G** A {u,J) = -Uj o dt J o due^-e^-i—^iu 

t 2 it 

— (qg s (? ■ Gq))[tpn{u) + 

16 q ■ v 



-((??) + tz(Q9s° ■ Gq))[vM + — g 2 {u) + t 2 g x (u)]} (21) 



+^f**J J I'duu J Vate'W-to^e^-^q.vfaiai) , (22) 

G^(u,u') = - 2 U r I f'due^-^e^lMu) + — <fe(u) + t 2 g x {u)]} 
it Jo t Jo q ■ v 

+ ^/t //« / ^e^^ (ai - UQ3), e l * (ai - m3) [v!|(a t ) - ^( a! ) + - «)^|(a,)] (23) 

where ^ = 1.76GeV, /„. = 132MeV, (qq) = -(225MeV) 3 , (qg s a ■ Gq) = m 2 (qq), m 2 = 
0.8GeV 2 . For large euclidean values of uj and u/ this integral is dominated by the region 
of small t, therefore it can be approximated by the first a few terms. 

After Wick rotations and making double Borel transformation with the variables u and 
u/ the single-pole terms in (|TTD are eliminated. Subtracting the continuum contribution 
which is modeled by the dispersion integral in the region u,u' > u c , we arrive at: 

9*Mk = - 5 J^e^ 4 {/i^ ff K)T 3 / 2 (f ) - a[^(u )Tf (f ) 

2 

-j>9i(uo) + tG 2 {uo)] + ^b w (tto) - ^gi(uo) + ^sG 2 (u )]} (24) 
+^/ 3 G WT 3 / 3 (f) , (25) 

n k 

where f n (x) = 1 — e~ x p- is the factor used to subtract the continuum, u c is the 

k=0 

continuum threshold. u = j^f-, T = ^+t 2 ' ^i' ^ 2 are ^ ne Borel parameters a = 
— (27r) 2 (gg). The functions G 2 {uq) and lf(uo) are defined as: 

G 2 (u ) = / g2{u)du , (26) 
Jo 
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jG, \ f ^ j r v?37r(«i, 1 + «o - 2ai,«i - m ) [ 1 ^("i, 1 - «i - o> 3 ,a 3 ) 

h ( u o) = / dai[ / da 3 = . 

Juo ai — u Jo a 3 

(27) 

We have used integration by parts to absorb the factors (q ■ v) and l/(q ■ v). In this 
way we arrive at the simple form after double Borel transformation. In obtaining (24) we 
have used the Borel transformation formula: B^e auJ = 5 (a — ^). 

Similarly we have: 

92fc*h = -^e^{^(u )T*f 3 (f)+4G 2 {u )T a f 1 (f) 

-4^( Mo )T 2 /i(f ) - f/^M + |gl/w(uo)} (28) 

+^ft G Wr 2 /i(|) , (29) 

where the function lf(u ) is defined as: 



lf(u ) = I " dat I —[<P\\(ati) - <p±(ai) + G~ — — —)<p\\(ati)] ■ (30) 



1 +"0 
2 



"11 



d« 3 _ 1 Qti — Mq, 



o a 3 ' 2 a 3 

From (2^) and ( P8| ) we know that both g 2 and (73 are negative using the notations in 
19] . In the following we will always discuss the absolute values of #2,3- 



3 Determination of the parameters 

In order to obtain the coupling constants from (p4T )-(P8l) we need the mass parameters A's 
and the coupling constants /'s of the corresponding interpolating currents as input. The 
results are |16| 

A A = 0.8 GeV f A = (0.018 ± 0.002) GeV 3 , 

A s = 1.0 GeV h = (0.04 ±0.004) GeV 3 . (31) 

For the sum rule fl2"4]) and (^3) the continuum threshold is uo c = (2.5 ± 0.1)GeV. 



We use the wave functions adopted in |12[ to compute the coupling constants. Moreover, 
we choose to work at the symmetric point T\ = T 2 = 2T, i.e., uq = | as traditionally done in 



literature ||12|| . Such a choice is very reasonable for the sum rules for g 2 since S* and S c are 
degenerate in the leading order of HQET. The mass difference between E* and A c is about 
0.2GeV. Due to the large values of 1\, T 2 ~ 3.2GeV 3> A used below, the choice of T\ = T 2 
is acceptable. We use the scale w = 1.3GeV, at which the values of the various functions 
appearing in (p4"D-(pS|) at Uq = ^ are: (p n (uo) = 1.22, ippiuo) = 1.142, f a (uo) = 1.463, 
gi (u ) = 0.034GeV 2 , G 2 (u ) = -0.02GeV 2 , / 3 g (m ) = -2.75 and lf(u ) = -0.24GeV 2 . 
We have used the asymptotic forms for the wave functions y? 37r (aj), f±(ai), <p\\(ai), <p±(o<i) 
and (p\\(oti) to calculate I 3 (u ) and lf(u ), since these wave functions are not known very 
well. f 37T = 0.0035GeV 2 . 

4 Numerical results and discussion 

We now turn to the numerical evaluation of the sum rules for the coupling constants. 
Since the spectral density of the sum rule (|24"D-([2"8|) p(s) is either proptional to s 2 or s 3 , 
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the continuum has to be subtracted carefully. We use the value of the continuum theshold 
u c determined from the corresponding mass sum rule at the leading order of a s and I/tuq 
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The lower limit of T is determined by the requirement that the terms of higher twists 
in the operator expansion is reasonably smaller than the leading twist, say < 1/3 of the 
latter. This leads to T > 1.3 GeV for the sum rules (p4"D -( p8|) . In fact the twist-four terms 
contribute only a few percent to the sum rules. The upper limit of T is constrained by 
the requirement that the continuum contribution is less than 50%. This corresponds to 
T < 2.2GeV. 

The variation of g 2 ^ with the Borel parameter T and u c is presented in Fig. 1 and 
Fig. 2. The curves correspond to uj c = 2.4, 2.5, 2.6GeV from bottom to top respectively. 
Stability develops for the sum rules and fl2"8p in the region 1.3 GeV <T<2.2 GeV, we 
get: 

92h*h = (2.5 ± 0.4) x 10~ 3 GeV 6 , (32) 
9sh*fA = (6.8 ± 0.4) x 10- 4 GeV 6 , (33) 

where the errors refers to the variations with T and u c in this region. And the central 
value corresponds to T = 1.6GeV and u c = 2.5GeV. 
Combining ([31]) we arrive at 

g 2 = 1.56 ±0.3 ±0.3 , (34) 
g 3 = 0.94 ±0.06 ±0.2 , (35) 

where the second error takes into account the uncertainty in f's. 

The recent CLEO measurement p0[ of the £* — > A c n decays gives g 3 = y / 3(0.57± 0.1) 



||2 1|| , where the factor \f?> arises from the different notaions. The decay £* — > S c 7r is 
kinematically forbidden so the direct measurement of g 2 is impossible. In the large-iV c 
limit of QCD (N c is the number of colors) g 2t3 is related to the nucleon axial charge 



gA = 1.25, g 2 = 3/2^, g 3 = y3/2<^4 |22|, p3|| . The quark model result is g 3 = y/3 x 0.61, 
#2 = 1-5 x (0.93 ± 0.16) I^TJ. The short-distance QCD sum rules with the external field 



method [fjl| yields g 2 = 1.5 x (0.4 ~ 0.7), g s = ^3/2 x (0.4 ~ 0.7). In this work we employ 
light-cone QCD sum rules to calculate the strong coupling constants #2,3- Both sum rules 
for g 2 and g 3 is very stable with reasonable variations of the Borel parameter T and the 
continuum threshold u c as can be seen from Fig. 1 and Fig. 2. It is interesting to note 
that the numerical values of g 2 and g 3 are consistent with both the experimental data and 



the quark model result ||21|| . It is also interesting to notice the deviation from the large 
N c limit prediction [22, p3fl . Moreover the result from the short-distance QCD sum rules 
is compatible with the present work with the light-cone QCD sum rule approach if we use 
the same values of u c and /'s though the errors are quite large. 
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Figure Captions 



Fig. 1. Dependence of gi on the Borel parameter T for different values of 
the continuum threshold tu c . From top to bottom the curves correspond 
to uj c = 2.6,2.5,2.4 GeV. 

Fig. 2. Dependence of g% on the Borel parameter T for different values of 
the continuum threshold u c . From top to bottom the curves correspond 
to lo c = 2.6,2.5,2.4 GeV. 
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